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Abstract 

We study the renormalization of the fermion mixing matrix in the Standard Model 
and derive the constraints that must be satisfied to respect gauge invariance to all 
orders. We demonstrate that the prescription based on the on-shell renormalization 
conditions is not consistent with the Ward-Takahashi Identities and leads to gauge de- 
pendent physical amplitudes. A simple scheme is proposed that satisfies all theoretical 
requirements and is very convenient for practical calculations. 



Despite a few interesting papers have been devoted to the fermion mixing renor- 
malization @, [|, |J, the subject has so far escaped much attention. This is mostly 
due to the fact that, as a result of GIM cancelations, the radiative corrections related 
to the renormalization of the CKM matrix can be made very small, 0(Gi?m^), where 
m q is the mass of a light quark. They are therefore of little practical importance in the 
context of the Standard Model (SM). Still, the subject has some conceptual interest 
in its own, not to mention the relevance of mixing in many extension of the SM. In 
this letter, we reconsider it from a different point of view, which allows us to point out 
some inconsistency in previous analyses and to propose an alternative solution. 

For definiteness, we concentrate on the case of the CKM matrix renormalization in 
the SM. A convenient framework to study this issue is provided by the Ward Takahashi 
Identities (WTI) of the theory with background fields @, ||, 0. Indeed, as the diag- 
onalization of the fermion mass matrix is achieved by field redefinitions that do not 
commute with the gauge transformations, the CKM elements appear explicitly in the 
WTI, unlike masses and gauge couplings. This will give us a strong constraint. At the 
functional level, the WTI which represent the a + generator of SU(2)l are implemented 
by the local functional operator W + acting on the effective action T (see M). For our 
purposes, the relevant part of W+ is the one which contains the quark fields: 



yyquark = 



)L yO _i §_ yO , L 

u ud 6$% 6i/>Z ud d 



(1) 



Here V° is the tree level CKM matrix. Upon renormalization, the fermionic fields 
are rescaled by non-diagonal complex wave function renormalization (WFR) matri- 
ces Z u and Zd\ consequently, V° is replaced by a renormalized CKM matrix V = 
(Z^)~^ V° (Z%)2, where Z L is the left-handed component of the WFR. Expanding V 
at first order, we obtain for the CKM counterterm 

SVua = \ E SZu*Vu> d ~ \ E VrtSZif (2) 

u' d! 

An additional constraint on the WFR comes from the requirement of unitarity for V, 
E V* (5Z% d + 8Z%) = E {SZuw + VW (3) 

d> u> 

Notice that this constraint follows also from the request that the commutation relations 
among the operators W+, W_, and W3 be preserved. Combining Eqs. (0) and (0) we 
find 

-7 TO(M«u-5Z^)V« i + y £ l V va (8ZS d -5Z£) , (4) 



ud 
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where, as expected for a unitary matrix, the renormalization of V is expressed in terms 
of anti-hermitian matrices. 

Any renormalization prescription that preserves the above WTI leads to a gauge- 
independent definition of the CKM matrix. To prove this theorem we start noticing 
that, as a consequence of W+(r) = 0, one has 



%vv + (r) = >v + (%r) + E 



a.d 



^ d ^wrwi dtV ^ l+ 



(5) 



where £ is a gauge parameter and the ellipses stand for additional contributions not 
relevant for us. On the other hand, the gauge variation of V is controlled by a Slavnov- 
Taylor Identity (STI), d^T = S(d x T) || ||, where x is the anti-commuting source of 
the composite operator generated by the variation of £, and S is the Slavnov- Taylor 
operator ||. As W+ does not depend on \ an d commutes with S, the first term of 
Eq. (H) vanishes. The only possibility compatible with the invariance of the theory 
(i.e. with Eq. (||])) is then that all the parameters in the square brackets of Eq. (f|) are 
identically zero. From this d^V = follows. 

An interesting point about Eq. @ is that the counterterms for physical parameters^, 
the CKM matrix elements, are given in terms of conventional objects like the WFR con- 
stants. From a rigorous point of view, the distinction between physical and conventional 
objects can be formulated in terms of the cohomology classes of the Slavnov- Taylor op- 
erator of the theory ||. The physical parameters on which the S- matrix depends are 
the coefficients of the cohomology classes. For example, in the absence of mixing, one 
cohomology class is provided by each Yukawa interaction term of the Lagrangian, which 
guarantees that the masses are physical objects. In the case of mixing, the Yukawa 
couplings are complex non-diagonal matrices and can be diagonalized by a redefinition 
of the fields, that is a finite WFR. In higher orders, however, the redefinition of fields 
originated by the diagonalization of the Yukawa matrix inevitably mixes with the one 
generated by the anomalous rescaling of the kinetic terms. A formal cohomological 
analysis [|7| shows that indeed the off-diagonal field redefinition contains some physical 
parameters, namely the CKM matrix elements. On the other hand, Eqs. allow us 
to disentangle the physical information related to the diagonalization of the Yukawa 
matrix (contained in the CKM matrix elements and in their counterterms) from the 
unphysical information carried by the Z factors. These constraints rely on the gauge 
invariance of the theory and, in a consistent framework, should be all satisfied. 

As we have seen, the definition of the CKM matrix at higher orders is conveniently 
expressed in terms of the anti-hermitian component of the WFR. It therefore depends 

1 Strictly speaking, the CKM element V u d is not a physical quantity; physical observables are 
identified at tree level by the Jarlskog reparametrization invariants fiofl , which are constructed from 
physical amplitudes. This scheme can in principle generalized at higher orders. 
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on the scheme chosen for the WFR. It should be clear, on the other hand, that once 
the counterterm 5V is calculated through Eq. @, it can be used independently of the 
choice of the SZ factors, because physical amplitudes are independent of the scheme 
adopted for the WFR. For example, in practical applications at the level of ^-matrix, 
it is often convenient to avoid the rescaling of the fields (i.e. the WFR) altogether 



TT | and introduce only the LSZ factors for the external fields: if mixing is present, 
however, one still has to renormalize the mixing parameters. 

A number of renormalization prescriptions for the CKM matrix are indeed possible; 
a first convenient option is the MS subtraction: by definition, assuming gauge invariant 
mass renormalization and after adjusting for the possible breaking of chiral invariance, 
it satisfies the WTI and the STL Hence, as a consequence of the above theorem, it can 



be guaranteed to yield a gauge-independent ultraviolet pole 5V-^ |12j to all orders (a 



proof can also be found following [13|). On the other hand, it is well-known that the 
decoupling of heavy particles is not manifest in the MS scheme. This means that if 
we work in the framework of an effective Lagrangian where the heavy fields (W boson 
and top quark) are integrated out, the dimension three and four operators that mix 
the quarks yield contributions to the amplitude which are not suppressed by the high 
mass scale. Moreover, as noted in |14| , these terms are not even defined in the limit of 



vanishing light quark masses. All this makes the MS definition unnatural in the context 
of effective Lagrangians. As the CKM elements are mostly determined from low-energy 
hadronic processes, this is not very convenient. Physical amplitudes calculated with an 



MS counterterm for the CKM would depend on the renormalization scale (see [|15| for 
studies of the scale evolution of the CKM matrix) and would contain 0(a) corrections 
proportional to (mf + m|)/(mf — mj), where rriij are the poorly known light quark 
masses. 

A second possibility consists in fixing four CKM elements in terms of four physical 
amplitudes, e.g. of the four most precise experimental processes. This procedure 
bypasses the definition of the WFR, but destroys the symmetry between the quark 
generations and is not practical in higher order calculations. 

A third option is provided by the use of the on-shell renormalization conditions of 



Ref-Hlql to define the WFR constants and, through Eq. (g), the CKM counterterm, as 
has been suggested by Denner and Sack || and generalized to extended models in j4|. 
This approach also implies decoupling in the sense explained before, i.e. dimension 
three and four operators are removed. The renormalization conditions, which are 
written in the u sector in terms of the fermionic two-point functions T m r(jf), 

u(m u ) T m r(f) = 0; T Uu/ (f) u(m u >) = 0; 

u(m u ) Tu U (jS)- = 1; Tuu(f)u(m u ) = 1, (6) 

i> - m u $ - m u 
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fix the masses of the u quarks and all the 5Z UU > by setting them equal to the LSZ factors. 
Their use in Eq. @ defines a counterterm SV u d which makes physical amplitudes finite 

We now show explicitly that the latter procedure leads to gauge dependent ampli- 
tudes in one-loop calculations. To this end, we consider the decay of a W boson into 
two arbitrary quarks u and d that was studied in Ref. ||. We conform to the notation 
of that paper and write the one-loop renormalized amplitude as 

M ud = V ud M (l + S vert + — - — + lsz w ] 
V e s w 2 J 

+M (\ £ 6ZX V U ' d + \ E V udl 5Z L d , d + 8V ud ) , (7) 

V u' d' / 

where A4o = —g/y/2 u(m u ) ^(M w ) a_ v (m^), u and v are the spinors of the final- 
state quarks, a± — (1 ± 7s)/2 are the right and left-handed projectors, and is the 
polarization vector of the W boson. 

Let us now consider the gauge dependence of the individual contributions to M. u d- 
For our purposes, it is sufficient to consider only the ^ w gauge parameter dependence. 
As the total amplitude in Eq. (|7|) must be gauge independent, d^ w M. u d = 0. Also the 



counterterms be and Ss w , as defined in [|TTJ, do not depend on £ w at the one-loop level. 



On the other hand, the gauge variation of the proper vertex 5 vert is not trivial; it can 
be studied using the STI that governs the gauge dependence of the Green functions^. 
On the mass shell, and after contracting with e M and with the external quark spinors, 
one finds the following non-linear identity written in terms of 1PI Green functions 

d( w T w + a d= E \^xn-w+^w+ud + ^XiUVu'^w+u'd + ^w+ud^XiV d 'd. • (8) 

i=l,2 



Here 7^ and r\ u ^ are the sources associated to the BRST variation of and u, d fields 
and xi,2 are the sources of the two independent composite operators generated by the 
variation of the two gauge fixing parameters £ Wj i and £w,2- Following the common 
practice, we have set £ w = = ^ w>2 - 
At the one-loop level, Eq. (W) reduces to 



u iw L W+ud ~ [K 2-^i 



r T(i) T/ v^rW t/ , v^t/ r (i) 



(9) 



where the superscript (1) indicates that the proper functions are evaluated in the one- 
loop approximation and T that only the transverse component is considered. All terms 
are evaluated on the mass-shell of the physical fields. In a similar way, one may find 

2 A detailed illustration of this kind of STI g, § will be given in ]TJ. 
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STI for the two-point functions of the W boson and of the quarks. At one-loop level 
and adopting the standard tadpole renormalization, which consists in removing the 
whole tadpole amplitude, one finds for the W boson WFR factor 







VT(1) 



9^ W SZ W — &£ w x. w+w . 



(P 2 



2Er 

i=l,2 



T(l) 



{Ml). 



(10) 



The treatment of the quark two-point functions is slightly more involved: in the case 
of the u quarks, for instance, we decompose the unrenormalized self-energy according 
to 

Sw(p) = s w $ a - + s w + s w +m u ,a + ). (11) 

The one-loop STI for £ uu / then reads 



9^E« (p) = - £ [r& v (jO (y - "v) + {i> - m u ) T% ]u M 



i=l,2 



(12) 



Splitting r^j^ into its left and right-handed components, r£-S^, we find on the 



L,R(1) 



mass shell of the u quark 

m u0^ w \^ uu i -f- ^ uu i 



p-R,(l) 



(13) 



One can then use the above equations in the definition of the on-shell 5Z L (see for ex. 
Eq. (3.4) of Ref.||) and find 



d^SZX = 2 £ r5«,(m u ). 



(14) 



i=l,2 



An analogous result holds in the d sector. 

Inserting Eqs. ( |T0l) and ( |T4D into Eq. (|9]), we can write the gauge dependence of the 
vertex as 



Vud ^vert 



-~ d^ w (v ud 5Z W + J2 8Z& V u , d + V ud , 5Z L d , d \ . (15) 

Z V u> d> J 



We observe that, as a consequence of the unitarity of V, in the last two terms of 
Eq. (|TJ) the dependence on the CKM matrix factorizes; this is consistent with the fact 
that the STI in Eq. (||) would be exactly the same, were the CKM matrix diagonal, and 
(from a diagrammatical point of view) that the one-loop vertices involve one and only 
one charged quark current. The limit for massless fermions of Eq. agrees with the 
expressions reported in |TS|] . Note also that all the gauge dependence of the vertex is 
contained in two-point function contributions. Using Eqs. (|7|) and (|T^), d^ w M. udi = 
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is reduced to d^ w 5V u d = 0. In other words, we have shown that the gauge-parameter 
dependence introduced by the LSZ external field factors (that coincide with the on- shell 
WFR constants) is completely absorbed by the proper vertex. On the other hand, V u d 
is a parameter of the bare Lagrangian and its renormalization condition, in the present 
framework, should preserve the gauge independence. 

It may therefore seem surprising that the counterterm defined on the basis of the 
on-shell conditions is gauge-dependent. Indeed, an explicit one-loop calculation yields 
for the £ w dependent part of the WFR in n dimensions 

6 u ' u ~^lV U ' d 'J (2nr(k*-Z w Ml)l(k + P y-ml,] ' (16) 

with p 2 = m\ and u 7^ u f ; again an analogous expression holds for SZ^ d ,. We see from 
Eq. (|l^) that the l/(n — 4) pole of the anti-hermitian part of 8Z L, ^ W is removed by 
GIM cancelations, ensuring that the divergence of 5V u d in Eq. (|j) is gauge independent. 
On the contrary, the momentum dependence of the integrand in Eq. (|16D spoils the 
GIM cancelations for the finite part of 5V u d in Eq. (|j), which turns then out to be 
gauge dependent. We conclude that the W-dec&y amplitude calculated in the on-shell 
framework of Ref. || is gauge dependent. This can be understood by noting that the 
finite part of the WFR factors defined on-shell violates the WTI, in particular it does 
not satisfy Eq. (|3]). 

A convenient and natural alternative to the prescription of ||, which maintains 
decoupling and enhances the symmetry among the quark generations, can be obtained 
by imposing the following conditions on the off-diagonal two-point functions: 



rw(o) = o ; rL'(o) ^^IWjQ a. 



0, (17) 



for u 7^ u', and analogously for the d sector. These conditions do not fix the diagonal 
two-point functions. They also do not fix the off-diagonal hermitian part of the right- 
handed WFR, which however is finite and suppressed by light quark masses. We 
choose to set SZ^f 1 = for u 7^ u' . The normalization point, $ = 0, is the same for all 
flavors and all divergences related to the mixing are correctly subtracted, as they are 
logarithmic and do not depend on masses and momenta. In addition, Eqs. (0) avoid 
problems in the treatment of the absorptive parts of the two-point functions whenever 
any of the quarks is heavy (typically, in the SM, this is the case with the top quark). 

At the one-loop level, the antihermitian and hermitian WFR factors 8Z U ^, and 
SZ^u' obtained from Eq. (|T7|) can be expressed in terms of self-energies evaluated at 
zero momentum transfer: 



SZ, 



L,A 



ml + ml, 



ml -ml, 



E^(0) + 2E^(0) ; **£f = -X&,(0). (1 
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It is straightforward to verify that SZ U ^, satisfies the condition of Eq. (§) because it 
reduces it to the case of no mixing. The CKM counterterm obtained by the use of 
5Z L ' A in Eq. (|) is, in units of g 2 / (64M^tt 2 ) , 



SVud — X! Vu'dVud'Vu'd' 
u'j^u d' 

n-5y 
2(y-i) ~ 



(ml + ml,) m\, 




e jj, 1 



d, vl «-> d') , (19) 



where y = m d ,/M^. and e = (4 — n)/2. SZ^, and 8V u d are gauge independent, as can 
be directly seen from Eq. (|T2|) . We stress once more that the use of this counterterm, 
based on Eq. ([[7]) and consequently on the Z L)A factors, is independent of the choice 
of WRF in the rest of a calculation, and corresponds to just one of the many gauge- 
invariant definition of the one-loop CKM elements. Needless to say, it is always the 
LSZ procedure (see Eqs. @) to dictate the treatment of the external lines. 5V u d can 
therefore be used without modifications in the calculation of ||: in that case, the 
results for the VT-decay amplitude are gauge independent and differ from || by gauge 
dependent 0(G^mf ight ) terms. 

We now consider the consistency of the renormalization conditions of Eqs. ( |T7|) 
beyond one-loop. First, we can show that the Eqs. (|17D respect the WTI for the Wud- 
vertex at all orders; this WTI reads 

iPu + Vd)n Fw+udiPu, Pd) ~ M w^G+ud(Pu, Pd) 

-^(j2 V ud>a + T s , d (p d ) ~Yl V u>d?uu>(-Pu)a^j = 0. (20) 

At the one-loop level and for on-shell amplitudes, Eq. ( |20| ) holds even in the case the 
external Goldstone boson G + and the W + are quantized. Differentiating with respect 
to p u and pa and setting all momenta to zero, one finds that 

5>«r r| d (0) - E r L(0) Vu'd = V ud G ad (0\ (21) 

d' u< 

where G a d(0) is a convergent term induced by Tg +iid . The second condition of Eqs. ( |17|) 
reduces this constraint to the case of no mixing at all orders. This is the crucial re- 
quirement. Similarly, the first condition of Eqs. (|I7|), used in Eq. ( pH|) at zero momenta, 
reduces it to the well-known constraint on the renormalization of the Yukawa coupling 
in the absence of mixing. 

We also investigate the effect of Eqs. (|T7| ) on the STI: at the one-loop level, they 
induce several constraints on the renormalization of x-dependent Green functions which 
appear in Eq. ([12]), e.g. rj^j ,; using them together with Eqs. ( |T7D in the two- loop STI, 
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the latter can be linearized and reduces to its one-loop form at p = 0. This is a non- 
trivial result, as the gauge dependence of the two-loop off-diagonal S„ M ,(p) can now 
be written only in terms of two-loop Green functions. At the one-loop level several 
simplifications occur: for example one has r^^,(0)a_ = 0, because only the left- 
handed source rj^, is involved. Beyond one-loop, r]$ may also contribute, to the effect 
that this and analogous simplifications are not granted any longer. Consequently, the 
combination S^ u ;(0) + 2 Sf u ;(0) is not guaranteed to be gauge invariant at two or more 
loops, although analogous but more complicated gauge-invariant combinations do exist. 
As already mentioned, what is certainly gauge invariant at all orders is the MS pole of 
5V; we have explicitly verified this property using the STL 

A related problem concerns E R , which is finite at one- loop just because of the 
GIM mechanism, but it may require a subtraction beyond one-loop; this would also 
modify the conditions of Eqs. (|1?D. Moreover, we note that a rigorous analysis of all the 
WTI (not only of Eq. (|20|)) at higher orders cannot be performed without specifying 
the whole set of normalization conditions. The previous points show that a complete 
discussion of a non-MS renormalization of the mixing matrix beyond one-loop becomes 
extremely complex [I7|, and is only partially simplified by Eqs. (|TTD . 

The renormalization conditions of Eqs. (|L7D can be used in any model containing 
Dirac fermion mixing. For example, all has been said applies directly to the case of 
lepton mixing, which is suggested by recent experiments, if the neutrinos have Dirac 
masses. They can also be generalized to extended models, along the lines of Refs. J3J. 

In summary, we have reanalyzed the renormalization of the fermion mixing param- 
eters in the Standard Model. We have reviewed several possibilities for the definition of 
the CKM matrix at higher orders, showing the constraints they have to satisfy in order 
to respect gauge invariance. In particular, we have demonstrated that the prescription 
based on the on- shell wave function renormalization constants is not consistent with the 
Ward-Takahashi Identities and leads to gauge dependent physical amplitudes. We have 
therefore proposed a simple scheme that naturally satisfies all theoretical requirements 
and is very convenient for practical calculations. 
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